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ABSTRACT 

A two-equation turbulence model of the K — s type was recently derived by Yakhot & 
Orszag based on Renormalization Group (RNG) methods. It was later reported that this 
RNG based model yields substantially better predictions than the standard K — e model for 
turbulent flow over a backward facing step - a standard test case used to benchmark the 
performance of turbulence models in separated flows. The apparent improvements obtained 
from the RN G K — e model were attributed to the better treatment of near wall turbulence 
effects. In contrast to these earlier claims, it is shown in this paper that the original version 
of the RNG K — e model substantially underpredicts the reattachment point in the backstep 
problem - a deficiency that is traced to the modeling of the production of dissipation term. 
However, with the most recent improvements in the RNG K — e model proposed by Yakhot 
and co-workers, excellent results for the backstep problem are now obtained. Interestingly 
enough, these results are not that sensitive to the details of the near wall treatment. 
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1. INTRODUCTION 

During the past two decades, the K — e model has played a central role in the calculation 
of many of the turbulent flows of scientific and engineering interest. In this two-equation 
model, the turbulence length and time scales are built up from the turbulent kinetic energy 
K and dissipation rate e which are obtained from separate modeled transport equations. The 
Reynolds stress tensor is represented by an eddy viscosity model constructed from length and 
time scales based on K and e. It was during the early 1970’s that the original version of the 
K — e model was developed by Launder and co-workers [1, 2] based on a phenomenological 
approach which incorporated empirically many of the earlier ideas of Kolmogorov [3] . In the 
intervening years, a variety of modifications have been made to the K — e model to account 
for other complicating features such as near- wall turbulence, body forces and compressibility 
effects to name a few [4-7]. 

One of the major criticisms that has been leveled against the K — e model is that it is 
ad hoc - namely, it is not derived from the Navier-Stokes equations in a systematic fash- 
ion. However, Yakhot and Orszag [8] recently derived a version of the K — e model using 
Renormalization Group (RNG) methods. In this approach, an expansion is made about an 
equilibrium state with known Gaussian statistics by making use of the correspondence prin- 
ciple wherein the effects of mean strains are represented by a random force. Bands of high 
wavenumbers (namely, the small scales) are systematically removed and space is rescaled 
in a manner analogous to that employed in. the study of phase transitions. The removal of 
only the smallest scales gives rise to subgrid scale models for large-eddy simulations whereas 
the removal of successively larger scales ultimately leads to Reynolds stress models. At high 
turbulence Reynolds numbers, the RNG based K — e model of Yakhot and Orszag [8] is of 
the same general form as the standard K — e model; however, the constants are calculated 
explicitly and assume somewhat different values. One major difference between the RNG 
and standard K — e model lies in the near wall treatment. The RNG model can be inte- 
grated directly to a solid boundary without the need for ad hoc wall damping functions. 
Several applications of this RNG K — e model have been subsequently reported by Yakhot, 
Orszag and co-workers. Probably the most notable one - since it deals with a complex shear 
flow with separation - is for turbulent flow over a backward facing step (see Karniadakis et 
al. [9]). It is well-known that the standard K — e model underpredicts the reattachment 
point in this problem - a deficiency that has been widely discussed in the literature since 
the 1980/81 Stanford Conference on Complex Turbulent Flows [10]. Karniadakis et al [9] 
claimed that the RNG K — e model predicts reattachment 7.3 step heights downstream of 
the. step corner for the Kim, Kline and Johnston [11] test case - a result that is extremely 
close to the experimental value. This apparent improvement in the predictions of the RNG 
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based K — e model were attributed largely to the better treatment of near wall turbulence 
effects [9]. 

Subsequent to the publication of the paper by Karniadakis et al. [9], several adjustments 
have been made in the RNG K — e model. An error in the calculation of the constant 
Cei in the dissipation rate transport equation was corrected (see Smith and Reynolds [12] 
and Yakhot and Smith [13]). Furthermore, a double expansion technique was introduced by 
Yakhot et al. [14] in an effort to yield improved models for the production of dissipation term 
that can better accommodate large strain rates. In the latter paper, the backstep problem 
for the Kim, Kline and Johnston [11] test case was re-computed yielding excellent results for 
this latest version of the RNG K — e model. Consequently, we feel that there is a need to 
clarify several questions: 


(1) Can the original version of the RNG K — e model yield good results for the backstep 
problem? 

(2) Does the latest version of the RNG K — e model yield excellent results for the backstep 
problem due to the change in the constant C s i, the addition of the new model for the 
production of dissipation term, or some other combination of factors? 

While the first question was answered to the negative in a recent paper by Thangam and 
Speziale [15], there are still some open questions regarding the near wall modeling that need 
to be addressed. In the sections to follow, an effort will be made to fully clarify these issues 
and to gain a better understanding of the role of the various modeled terms in the prediction 
of turbulent separated flows. 

2. THE RNG K-e MODEL AND THE BACKSTEP PROBLEM 


The RNG procedure of Yakhot and Orszag [8] for incompressible turbulent flows yields 
the renormalized equations of motion 


dui _ dui __ 
dt dxj 


dP_ 

dxi 


d_ ' 

dxj _ 


(u -t- u T ) 




dui 

dxi 


= 0 


( 1 ) 

( 2 ) 


where «,• is the mean velocity, P is the mean pressure, v is the molecular viscosity and vj is 
the eddy viscosity. At high turbulence Reynolds numbers, the eddy viscosity takes the form 


^ I<2 
VT = C, T 


( 3 ) 
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where 


(4) 


1 - 




K = — u'-ui-, 
2 1 •’ 


e = v 


du'i du\ 
dxj dxj 


are, respectively, the turbulent kinetic energy and dissipation rate. is a dimensionless 
constant which was computed to be 0.085 using the RNG methodology. The turbulent 
kinetic energy and dissipation rate are typically obtained from modeled transport equations 
which, at high Reynolds numbers, take the form 


m_ dK- _d_ f^dK_' 

dt Ul dxi £ dxi l <tk dxi 
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where V is the turbulence production given by 


and 


V = 2uxSijSij 


^ 1 (dui . dm 
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( 8 ) 


vj dxij 

is the mean rate of strain tensor. The modeled transport equations (5) - (6) contain four 
coefficients: C e i, C s 2 , cr % and a e . In the standard K — e model, these coefficients are con- 
stants which are obtained from benchmark experiments for equilibrium turbulent boundary 
layers and isotropic turbulence. On the other hand, in the original version of the RNG K — e 
model of Yakhot and Orszag [8], the coefficients are constants which are calculated explicitly 
by the theory. In the newest version of the RNG K — e model developed by Yakhot et al, 
[14], corrections to the constants in the dissipation rate equation made by Yakhot and Smith 
[13] were implemented along with a modification of the production of dissipation term to 
account for large strain rates. These coefficients can be summarized as follows: 


Standard K — e Model [ 2 ] 


C„ = 0.09, C e 1 = 1.44, Ce 2 = 1.92 
&K = 1.0, (T e = 1.3 

RNG K-e Model [ 8 ] 

C „ = 0.085, Cei = 1-063, Cei = 1.72 
o K - 0.7179, o e = 0.7179 


(9) 


( 10 ) 
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New RNG K-e Model [14] 


C„ = 0.085, C A = 1.42 - 

C e 2 = 1.68, a K = 0.7179, <r £ = 0.7179 


(H) 


where rj = SK/e , S = (2S ij S ij )^ 2 , tj 0 = 4.38 and j3 = 0.015. The new RNG K-e 
model is obtained by implementing the corrections of Yakhot and Smith [13] (wherein C e \ 
was changed from 1.063 to 1.42 and C e 2 from 1.72 to 1.68) along with a strain-dependent 
modification to C e i which is significant for flows with large strain rates. In so far as the 
latter modification is concerned, y 0 is the fixed point for homogeneously strained turbulent 
flows and f} is a constant which was evaluated to yield a von Karman constant k ph 0.4 (see 
Yakhot et al. [14]). 

The test case to be considered is turbulent flow over a backward facing step (see Figure 
1). Model predictions will be compared with the experimental data of Kim, Kline and 
Johnston [11] as updated by Eaton and Johnston [16] for the 1980/81 Stanford Conference 
on Complex Turbulent Flows. For this flow configuration, the expansion ratio E (step 
height: outlet channel height) is 1:3 and the Reynolds number Re is 132,000 based on the 
inlet centerline mean velocity and the outlet channel height. 

The fully- developed mean velocity is of the two-dimensional form 

u = H(x,y)i + v(x,y)j. (12) 


(13) 


(14) 


Hence, the governing field equations to be solved are given by 
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where the turbulence production V takes the form 
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given that the eddy viscosity v? is as defined in (3). 

Equations (13) - (17) are solved subject to the following boundary conditions: 


( 16 ) 


(17) 


(18) 


(a) The inlet mean velocity profile u is specified five step heights upstream of the step 
corner. This profile is matched to the experimental data [16] by performing a separate 
developing channel flow calculation which also provides companion profiles for K and 
e at the inlet. 


(b) Extrapolated outflow conditions are applied thirty step heights downstream of the step 
corner. 


ti + = 


y + , 


(c) The law of the wall is applied at the upper and lower walls in the three layer form [15]: 

for y + < 5 

-3.05 + 5 in y + , for 5 < y + < 30 (19) 

[ 5.5 + 2.5 in y + , for y + > 30 


where y + = yu r / v, u + = uju T and u T is the friction velocity constructed from the wall 
shear stress in the standard manner. The law of the wall (18) is used in conjunction 
with the boundary condition of vanishing normal derivative of K at the wall along 
with the log layer formulas 


K_ 

u 2 t 



e = cy 


K-l‘ 

ny 


( 20 ) 


(where k « 0.4 is the von Karman constant) which are suitably interpolated to establish 
boundary conditions for K and e at the first grid point away from the wall [17]. The 
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friction velocity u r is approximated by its log law form u T « C X J A K X !' 1 in order to 
avoid problems with singularities at the reattachment point where the wall shear stress 
vanishes. 

The equations of motion are solved in their time dependent form by a finite volume 
method. A variable 200 x 100 mesh, which concentrates mesh points near the step corner, 
is used. This resolution is required to ensure the numerical accuracy of the results (many 
earlier computations of the backstep problem were under-resolved; see Thangam and Speziale 

[15] for a discussion of this issue as well as for other details of the numerical method). The 
computations were conducted on a CRAY-YMP computer. Approximately 15 minutes of 
CPU time is required to obtain a steady state solution for a given test case. 

3. DISCUSSION OF THE RESULTS 

First, we will present computed results for the standard K — e model since this is an im- 
portant benchmark against which the other results will be judged. All model predictions will 
be compared with experimental data for the Kim, Kline and Johnston test case as updated by 
Eaton and Johnston [16] for the 1980/81 Stanford Conference on Complex Turbulent Flows. 
In Figures 2(a)-(b) the mean velocity streamlines and mean velocity profiles obtained from 
the standard K — e model are compared with the experimental data. The standard K — e 
model predicts reattachment at Xr/H ph 6.25 - a value that constitutes approximately a 
12% underprediction of the mean experimental reattachment point of Xr/H « 7.1. Earlier 
reported results [10] of a 20-25% underprediction of the reattachment point appear now to 
have been in error due to insufficient numerical resolution as discussed by Thangam and 
Speziale [15] and Avva, Kline and Ferziger [18]. 

One of the criticisms that has been leveled against most of the K — e model calculations 
for the backstep problem lies in the use of wall functions. It is well-known that wall functions 
do not formally apply to separated flows. However, since the separation point is fixed - and, 
perhaps, since a large proportion of the turbulent kinetic energy is not associated with the 
separated zone - it now appears that wall functions do not give rise to major errors in 
the backstep problem, unlike in other separated flows. A systematic study of near wall 
turbulence models by the authors has indicated that they generally do not give rise to more 
than a 5% improvement in the results for the backstep problem obtained from well-resolved 
computations using wall functions. In Table 1, the reattachment point prediction of the 
standard K — e model using three-layer wall functions is compared with experimental data 

[16] as well as with the result obtained using the asymptotically consistent near wall model 
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of Speziale, Abid and Anderson [19]. Here, this particular near wall model is chosen since it 
was shown by those authors to be better behaved than four independent near- wall turbulence 
models that were recently reviewed by Patel, Rodi and Scheuerer [20] for boundary layer 
flows. It is clear from Table 1 that the use of an asymptotically consistent near wall model 
does not lead to a significant improvement in the results obtained using wall functions. 
Furthermore, it should be noted that the use of an asymptotically inconsistent near wall 
model (some of which are discussed in Patel, Rodi and Scheuerer [20]) can lead to even 
worse results than those obtained using wall functions! Consequently, we believe that it is 
acceptable to use wall functions when comparing the predictions of a variety of turbulence 
models in the backstep problem. 

Now we will present computed results for the original version of the RNG K — e model 
(Yakhot and Orszag [8]). In Figures 3(a) - (b) the computed mean velocity streamlines 
and mean velocity profiles are compared with the experimental data of Eaton and Johnston 
[16]. The surprising finding is that the predicted mean reattachment point of Xr/ I i ~ 4.0 
constitutes over a 40% underprediction of the mean experimental result of Xr/ H « 7.1. 
This deficiency — as well as a plausible explanation for its cause — was pointed out recently 
by Thangam and Speziale [15]. In a homogeneous shear flow, the eddy viscosity grows 
exponentially: 

Vt ~ exp(Ai*) 


where t* is the time (nondimensionalized by the shear rate) and A is the growth rate given 
by (see Speziale [21]) 


C„(Cg ~ i ) 2 

l(c„ - 1)(C. 2 - 1)J 


1/2 


( 21 ) 


which becomes singular for C e \ = 1. Consequently, for C e \ = 1.063, the growth rate of the 
eddy viscosity will be overly large. When a model overpredicts the eddy viscosity it will 
be too dissipative - a deficiency that will lead to an underprediction of the separated flow 
region in the backstep problem. Earlier reported results by Karniadakis et al. [9] for the 
backstep problem - which suggested that the model predictions were excellent - now appear 
to have been in error. 

With the correction of C eX to 1.42 by Yakhot and Smith [13], the problem of overpredicting 
the growth rate of homogeneous shear flow is eliminated. Speziale, Gatski and Fitzmaurice 
[22] showed that with the new value of C eX , the RNG K — e model yields excellent results 
for homogeneous shear flow. However, an additional production term in the dissipation 
rate transport equation was uncovered which Yakhot and Smith [13] were not able to close 
systematically. The new term was shown to be important for large strain rates. If this 
term is neglected, a revised RNG I< — e model is obtained which is of the same general 
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form as the Yakhot and Orszag [8] model; only the values of the constants are altered. In 
Figures 4(a)-(b), the mean velocity streamlines and mean velocity profiles obtained from the 
Yakhot-Smith RNG K — e model (ignoring the additional production term) are shown. It 
is clear that the model drastically overpredicts the reattachment point (i.e., Xr/H « 9.7 in 
comparison to the mean experimental result of 7.1). The origin of this erroneous prediction 
is clear: with the additional production term in the e-transport equation neglected, the 
Yakhot-Smith model yields a von Karman constant of 0.23 instead of the traditional value 
of 0.4 (see Yakhot and Smith [13]). This strongly indicates that the additional production 
term cannot be neglected in wall bounded turbulent flows. 

Yakhot et al. [14] recently developed a model for the additional production term in the 
e-transport equation by means of a scale expansion. More precisely, an expansion of this 
term was made in the ratio of the turbulent to mean strain time scale r/ = SK/e , where 
S = (2SijSijY^ 2 is the norm of the mean rate of strain tensor. The interesting finding was 
that no finite truncation of the expansion for the production of dissipation term in powers of 
rj suffices - terms of all orders must be retained to satisfy the crucial weak and strong strain 
limits [14]. This complication eliminated the possibility of obtaining a closed form solution 
for this term. A Pade approximation was made which had one undetermined constant (3 that 
was evaluated empirically by setting the von Karman constant to 0.4. This leads to a net 
production of dissipation term of the traditional form C s x(e/K)V with a variable coefficient 
given by 

ft. = 1-42 - (22) 

where f3 = 0.015 and r/o = 4.38 is the fixed point for equilibrium homogeneous flows. The 
addition of the second term to the r.h.s. of (22) - which becomes significant for large strain 
rates rj > r) 0 - is the feature that distinguishes the Yakhot et al. [14] model from the Yakhot 
and Smith [13] model. 

In Figures 5(a)-(b), the computed mean velocity streamlines and mean velocity profiles 
for the Yakhot et al. [14] RNG K — e model are compared with experimental data [16]. The 
model predicts reattachment at Xr/H & 6.7: a result that is within 5% of the experimental 
mean reattachment point of 7.1. Most notably, this new RNG K — e model - unlike all of 
the other models - yields a significant secondary separation bubble below the corner of the 
step that is more in line with the experimental data [23]. Yakhot et al. [14] showed that 
when this model is extended to include an anisotropic eddy viscosity, the predicted mean 
reattachment point is almost identical to the experimental value of Xr/H ^ 7.1. 
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4. CONCLUDING REMARKS 


Three independent versions of the RNG K — e model have been applied to the test case 
of turbulent flow over a backward facing step in order to assess its performance in separated 
flows. The main objective of this study was to better understand the impact that the recent 
changes in the Yakhot k Orszag RNG K — e model have had on the predictive capabilities of 
the model - particularly in a complex shear flow of engineering interest. It was found that the 
original version of the RNG K — e model substantially underpredicts the reattachment point 
for the Kim, Kline and Johnston test case due to an overprediction of the eddy viscosity 
(the model is far too dissipative). However, the latest version of the RNG K — e model 
[14] - which is an extension of the Yakhot and Smith [13] model to accommodate large 
strain rates in the production of dissipation term - yields excellent results for the backstep 
problem. The mean reattachment point is predicted to within 5% of the experimental value 
in comparison to the standard K — e model which yields a 12% error. Furthermore, a 
secondary separation bubble was obtained that is much more in agreement with that which 
is observed in experiments. These results are strongly dependent on the presence of the 
variable strain rate term in (21) and the calculations vividly demonstrate the sensitivity of 
the reattachment point to the values of C s \ and C e i. The encouraging point, however, is that 
the results are not that sensitive to the details of the near wall treatment and are obtained 
with no ad hoc adjustments of the constants. Based on these results, we feel that the current 
version of the RNG K—e model can be a useful turbulence model for practical engineering 
and scientific calculations. 
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NEAR- WALL MODEL 

Xr/H 

Three-layer wall functions 

6.25 

Speziale et al. [19] 

6.40 

Experimental data [16] 

7.1 


Table 1. The reattachment point in the backstep problem: Comparison of the predictions of 
different near wall corrections to the standard K — e model with experimental data. 
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X 


Figure 1. Schematic of turbulent flow over a backward facing step. 


12 


— K-e Model 

o Experimental Data [16] 



(a) Streamlines 



(b) Dimensionless Mean Velocity Profiles 


Figure 2. Comparison of the predictions of the standard K — e model [2] with experimental 
data [16]. (a) Streamlines and (b) mean velocity profiles. 
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— K-e Model 


o Experimental Data [16] 



(a) Streamlines 



(b) Dimensionless Mean Velocity Profiles 


Figure 3. Comparison of the predictions of the RNG K — e model of Yakhot and Orszag [8] 
with experimental data [16]. (a) Streamlines and (b) mean velocity profiles. 


14 



K-e Model 


o Experimental Data [16] 



(a) Streamlines 


X/H : 1.33 2.67 5.33 6.22 8.00 10.67 



U 

4 > 

(b) Dimensionless Mean Velocity Profiles 


Figure 4. Comparison of the predictions of the RNG K — e model of Yakhot and Smith [13] 
- where the additional production of dissipation term is neglected - with experimental data 
[16]. (a) Streamlines and (b) mean velocity profiles. 
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mo eling of the production of dissipation term. However, with the most recent im- 
provements in the RNG K-e model proposed by Yakhot and co-workers, excellent results 
or the backstep problem are now obtained. Interestingly enough, these results are 
not that sensitive to the details of the near wall treatment. 
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